Antibunched photons emitted by a quantum point contact out of equilibrium 
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Motivated by the experimental search for "GHz nonclassical light", we identify the conditions 
under which current fluctuations in a narrow constriction generate sub-Poissonian radiation. An- 
tibunched electrons generically produce bunched photons, because the same photon mode can be 
populated by electrons decaying independently from a range of initial energies. Photon antibunching 
becomes possible at frequencies close to the applied voltage V x e/h, when the initial energy range 
of a decaying electron is restricted. The condition for photon antibunching in a narrow frequency 
interval below eV/h reads E„T"n(l ~ T ™)T < 2^„[T„(1 - T n )] 2 , with T„ an eigenvalue of the 
transmission matrix. This condition is satisfied in a quantum point contact, where only a single T n 
differs from or 1. The photon statistics is then a superposition of binomial distributions. 
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In a recent experiment pj, Gabelli et al. have mea- 
sured the deviation from Poisson statistics of photons 
emitted by a resistor in equilibrium at mK tempera- 
tures. By cross-correlating the power fluctuations they 
detected photon bunching, meaning that the variance 
Varn = (n 2 ) — (n) 2 in the number of detected photons 
exceeds the mean photon count (n) . Their experiment is 
a variation on the quantum optics experiment of Hanbury 
Brown and Twiss [2|, but now at GHz frequencies. 

In the discussion of the implications of their novel ex- 
perimental technique, Gabelli et al. noticed that a general 
theory |3| for the radiation produced by a conductor out 
of equilibrium implies that the deviation from Poisson 
statistics can go either way: Super-Poissonian fluctua- 
tions (Var n > (n) , signaling bunching) are the rule in 
conductors with a large number of scattering channels, 
while sub-Poissonian fluctuations (Varn < (n), signal- 
ing antibunching) become possible in few-channel con- 
ductors. They concluded that a quantum point contact 
could therefore produce GHz nonclassical light 

It is the purpose of this work to identify the conditions 
under which electronic shot noise in a quantum point 
contact can generate antibunched photons. The physical 
picture that emerges differs in one essential aspect from 
electron-hole recombination in a quantum dot or quan- 
tum well, which is a familiar source of sub-Poissonian 
radiation 0, IE ■ ^ n those systems the radiation is pro- 
duced by transitions between a few discrete levels. In a 
quantum point contact the transitions cover a continuous 
range of energies in the Fermi sea. As we will see, this 
continuous spectrum generically prevents antibunching, 
except at frequencies close to the applied voltage. 

Before presenting a quantitative analysis, we first dis- 
cuss the mechanism in physical terms. As depicted in 
Fig. ^ electrons are injected through a constriction in 
an energy range eV above the Fermi energy Ep, leaving 
behind holes at the same energy. The statistics of the 
charge Q transferred in a time t h/eV is binomial 
with VarQ/e < (Q/e). This electron antibunching is 
a result of the Pauli principle. Each scattering channel 
n = 1, 2, . . .N in the constriction and each energy inter- 



val SE — h/r contributes independently to the charge 
statistics. The photons excited by the electrons would 
inherit the antibunching if there would be a one-to-one 
correspondence between the transfer of an electron and 
the population of a photon mode. Generically, this is not 
what happens: A photon of frequency oj can be excited 
by each scattering channel and by a range eV — Huj of 
initial energies. The resulting statistics of photocounts is 
negative binomial Q , with Var n > (n) . This is the same 
photon bunching as in black-body radiation [9(. 

In order to convert antibunched electrons into an- 
tibunched photons, it is sufficient to ensure a one-to- 
one correspondence between electron modes and photon 
modes. This can be realized by concentrating the current 
fluctuations in a single scattering channel and by restrict- 
ing the energy range eV — huj. Indeed, in a single-channel 
conductor and in a narrow frequency range lo < eV/% we 
obtain sub-Poissonian photon statistics regardless of the 
value of the transmission probability. In the more gen- 
eral multi-channel case, photon antibunching is found if 
[£„T„(1-T n )] 2 < 2J2 n [T n (l~T n )} 2 (withT n aneigen- 
value of the transmission matrix product tV). 

Starting point of our quantitative analysis is the gen- 
eral relationship of Ref. [3] between the photocount dis- 
tribution P(n) and the expectation value of an ordered 
exponential of the electrical current operator: 
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We summarize the notation. The function F(^) = 
Sfclo(£ fc /^')( nfc }f ^ s t ne generating function of the fac- 
torial moments (n k )f = (n(n — l)(n — 2) • • • (n — k + 1)). 
The current operator I = I out — L m is the difference of 
the outgoing current J out (away from the constriction) 
and the incoming current L ln (toward the constriction). 
The symbol O indicates ordering of the current operators 
from left to right in the order V m , /| ut , 7 out , L m . The real 
frequency-dependent response function "f(uj) is propor- 
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FIG. 1: Schematic diagram of a constriction in a conduc- 
tor (bottom) and the energy range of electronic states (top), 
showing excitations of electrons (black dots) and holes (white 
dots) in the Fermi sea. A voltage V drops over the con- 
striction. Electrons (holes) in an energy range eV — hcj can 
populate a photon mode of frequency ui, by decaying to an 
empty (filled) state closer to the Fermi level. 



tional to the coupling strength of conductor and photodc- 
tector and proportional to the detector efficiency. Posi- 
tive (negative) lo corresponds to absorption (emission) of 
a photon by the detector. We consider photodetection by 
absorption, hence 7(w) = for to < 0. Integrals over fre- 
quency should be interpreted as sums over discrete modes 
Lu p = p x 27t/t, p = 1,2,3,.... The detection time r is 
sent to infinity at the end of the calculation. We denote 
7 p = 7(w p ) x 2tt/t, so that J du>j(uj) ~ > J2 P lp- For ease 
of notation we set % = 1, e = 1. 

The exponent in Eq. is quadratic in the current 
operators, which complicates the calculation of the ex- 
pectation value. We remove this complication by intro- 
ducing a Gaussian field z{lj) and performing a Hubbard- 
Stratonovich transformation, 
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The brackets (■ ■ •) now indicate both a quantum mechani- 
cal expectation value of the current operators and a clas- 
sical average over independent complex Gaussian vari- 
ables z p = z(ujp) with zero mean and variance (|z p | 2 ) = 

1/7p- 

We assume zero temperature, so that the incoming cur- 
rent is noiseless. We may then replace I by J out and re- 
strict ourselves to energies e in the range (0, V) above 
E F . Let b] n (e) be the operator that creates an outgoing 
electron in scattering channel n at energy s. The out- 
going current is given in terms of the electron operators 
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deY,bi(e)b n (e + uj). 
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Energy e p — p x 27r/r is discretized in the same way as 
frequency. The energy and channel indices p, n are col- 
lected in a vector b with elements b pn — (27r/r) 1 / 2 6 n (e p ). 
Substitution of Eq. © into Eq. © gives 
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The exponents contain the product of the vectors 
b, w and a matrix Z with elements Z pn y n i = 
^nn' z p-p'lp-p' ■ Notice that Z is diagonal in the chan- 
nel indices n, n' and lower-triangular in the energy indices 
P,P'- 

Because of the ordering O of the current operators, the 
single exponential of Eq. JSJ factorizes into the two non- 
commuting exponentials of Eq. JSJ. In order to evaluate 
the expectation value efficiently, we would like to bring 
this back to a single exponential — but now with normal 
ordering J\f of the fcrmion creation and annihilation op- 
erators. (Normal ordering means b^ to the left of b, with 
a minus sign for each permutation.) This is accomplished 
by means of the operator identity |lfjj | 
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valid for any set of matrices Ai . The quantum mechanical 
expectation value of a normally ordered exponential is a 
determinant [Tl|. 
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Det(l + AB), B i:i = (b]b i ). (7) 



In our case A = e z e z - 1 and B = tt\ with t the TV x N 
transmission matrix of the constriction. 

In the experimentally relevant case Q, 0] the response 
function 7(0;) is sharply peaked at a frequency fl < V, 
with a width A <C ri. We assume that the energy depen- 
dence of the transmission matrix may be disregarded on 
the scale of A, so that we may choose an e-independent 
basis in which ttf is diagonal. The diagonal elements 
are the transmission eigenvalues Ti, T2, . . . Tjv <E (0,1). 
Combining Eqs. (JBUZI we arrive at 
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(In the second equality we used that Det e z+ = 1 , since 
Z is a lower-triangular matrix.) The remaining average 
is over the Gaussian variables z p contained in the matrix 
Z. 
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Since the interesting new physics occurs when f2 is close 
to V, we simplify the analysis by assuming that 7 (w) = 
for uj < V/2. For such a response function one has Z 2 = 
0. (This amounts to the statement that no electron with 
excitation energy e < V can produce more than a single 
photon of frequency uj > V/2.) We may therefore replace 
e z -> 1 + Z and e~ zt -> 1 in Eq. ©. We then apply 
the matrix identity 

Det (1 + A + B) = Det (1 — AB), if A 2 = = B 2 , (9) 

and obtain 
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We have defined (,X = ZZ^ and written out the Gaussian 
average. The Hermitian matrix X has elements 
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The integers p, p' , q range from 1 to Vt/2tt. 

The Gaussian average is easy if the dimensionless shot 
noise power S = ^ n T n (l — T n ) is 3> 1. We may then do 
the integrals of Eq. i|10|) in saddle-point approximation, 
with the result (13 

lnF(0 = ~ £ duj In [1 - £S 7 (w)(V - w)] . (12) 

The logarithm lnF(£) is t he g enerating function of the 
factorial cumulants ((n k }}{ [14(. By expanding Eq. I|12|) 
in powers of £ we find 
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Eqs. I|12ll and (|13|l represent the multi-mode superpo- 
sition of independent negative-binomial distributions [9j . 
All factorial cumulants are positive, in particular the sec- 
ond, so Varrt > (n). This is super-Poissonian radiation. 

When S is not 3> 1, e.g. when only a single chan- 
nel contributes to the shot noise, the result ijl2ll3(l re- 
mains valid if V — Q 3> A. This was the conclusion of 
Ref. 0, that narrow-band detection leads generically to 
a negative-binomial distribution. However, the saddle- 
point approximation breaks down when the detection 
frequency f2 approaches the applied voltage V. For 
V — SI < A one has to calculate the integrals in Eq. 
iJTUj l exactly. 

We have evaluated the generating function (|10[1 for a 
response function of the block form 



7o if V - A < uj < V, 
if uj < V - A, 
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with A < V/2. The frequency dependence for uj > V is 
irrelevant. In the case N = 1 of a single channel, with 
transmission probability T\ = T, we find 0] 

= — / duj\n[l + ^T(l-T)(V-uj)} 
tA (1 + x) ln(l +x)-x 
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with x = ^7oT(l — T)A. This is a superposition of bino- 
mial distributions. The factorial cumulants are 
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The second factorial cumulant is negative, so Varn < 
(ri). This is sub-Poissonian radiation. 

We have not found such a simple closed-form expres- 
sion in the more general multi-channel case, but it is 
straightforward to evaluate the low-order factorial cumu- 
lants from Eq. (jTUjl . We find 
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»f = ^(7oA) 3 - (35? - 155x5 2 + 155 3 ) ,(19) 
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with S p = $Z_[2n(l — T n )] p . Antibunching therefore re- 
quires S 2 < 2S 2 . 

The condition on antibunching can be generalized to 
arbitrary frequency dependence of the response function 
-f(uj) in the range V — A < uj < V of detected frequencies. 
For A < V/2 we find 
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x (2S 2 - AS 2 -{V- u)S 2 -^j 7 (w). (20) 

We see that the antibunching condition S 2 < 2S 2 de- 
rived for the special case of the block function (|14|) is 
more generally a sufficient condition for antibunching to 
occur, provided that dj/duj > in the detection range. 
It does not matter if the response function drops off at 
uj > V, provided that it increases monotonically in the 
range (V — A, V). A steeply increasing response function 
in this range is more favorable, but not by much. For 
example, the power law 7(0;) oc (uj — V + A) p gives the 
antibunching condition S 2 < 2S 2 x [1 + p/ (1 +p)], which 
is only weakly dependent on the power p. 

In conclusion, we have presented both a qualitative 
physical picture and a quantitative analysis for the con- 
version of electron to photon antibunching. A simple 
criterion, Eq. I|18(l . is obtained for sub-Poissonian pho- 
ton statistics, in terms of the transmission eigenvalues 
T n of the conductor. Since an iV-channel quantum point 
contact has only a single Tjv different from or 1, it 
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should generate antibunched photons in a frequency band 
(V — A, V) — regardless of the value of T N . The statis- 
tics of these photons is the superposition l|15(l of binomial 
distributions, inherited from the electronic binomial dis- 
tribution. There are no stringent conditions on the band 
width A, as long as it is < V/2 (in order to prevent multi- 
photon excitations by a single electron ^(|)- This should 



make it feasible to use the cross-correlation technique of 
Ref. to detect the emission of nonclassical microwaves 
by a quantum point contact. 
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